Mon. Not. R. Astron. Soc. OOP. [THSTI (20111 Printed 27 July 2011 (MN I^TeX style file v2.2) 



Are Halo and Galaxy Formation Histories Correlated? 



Jeremy L. Tinker^, Andrew R. Wetzel^, and Charlie Conroy^ 

^ Center for Cosmology and Particle Physics, Department of Physics, New York University, New York, NY 
^ Department of Astronomy, Yale University, New Haven, CT 
■^ Harvard-Smithsonian Center for Astrophysics, Combridge, MA 



o 

(N 



in 

(N 

u 

6 

^\ 



> 

o 
vn 

i> 

o 



% 



27 July 2011 



ABSTRACT 

The properties of dark matter halos, including mass growth, correlate with larger 
scale environment at fixed mass, an effect known as assembly bias. However, whether 
this environmental dependence manifests itself in galaxy properties remains unclear. 
We apply a group-finding algorithm to Data Release 7 of the Sloan Digital Sky Survey 
to estimate the dark matter halo mass of each galaxy and to decompose galaxies 
into those that exist at the centers of distinct halos and those that orbit as satellites 
within larger halos. Using the 4000-A break as a measure of star formation history, we 
examine the correlation between the quenched fraction of galaxies, /q, and large-scale 
environment, p. At all galaxy magnitudes, there is a positive, monotonic relationship 
between /q and p. We use the group catalog to decompose this correlation into the 
contribution from central and satellite galaxies as a function of halo mass. Because 
satellites are more likely to be quenched than central galaxies, the observed /q-p 
correlation is primarily due to variations of the halo mass function with environment, 
which causes a larger fraction of satellite galaxies at high p. For high-mass central 
galaxies, there is a weak but positive correlation between /q and p at both fixed galaxy 
luminosity and fixed halo mass. For low-mass central galaxies (M* < IO^^^Mq/Zi^), 
there is no correlation between /q and p. The latter results are inconsistent with the 
strong assembly bias of dark matter halos seen in this mass regime if recent galaxy 
growth at all correlates with recent halo growth, as we demonstrate through a high- 
resolution A^-body simulation. We also find that the mean stellar age of quenched 
central galaxies is independent of p at fixed M* , while the formation times of low mass 
halos vary significantly. We conclude that the processes that halt the star formation 
of low mass central galaxies are not correlated to the formation histories of their host 
halos, and old galaxies do not reside preferentially in old halos. 

Key words: cosmology: observations — galaxies: clustering — galaxies: groups: general 
— galaxies: clusters: general — galaxies: evolution 



1 INTRODUCTION 

Many galaxy properties are correlated with their large- 
and small-scale environment. The primary examples are 
the color-density rela tion a nd the m orphology-density re- 
lation (e.g., [Ocmlcr 1974 iDavis fc G cUcr 1976; Dr cssleij 
1 198(3 ). These results demonstrated that, in the local Uni- 
verse, galaxies in lower-density regions are typically less 
massive, and that even at fixed luminosity or mass, they 
are bluer and more likely to be spirals. With the advent of 
large-scale galaxy redshif t surveys such as the Sloan Dig- 
ital Sky Survey (SDSS; lYork et all I2OO0I '). these seminal 
works have been r e fined with increased stat i stical precision 
jHogget all |2004 



I2OO5I : iBaldrv et all 



Kauffmann et al 
20061 : 



Park et al 



I2OO4I : iBlanton et al 
I2OO7I : iBamford et~ 



[2OO9I) and extended via direct measurements of the cor- 
relations between s tar formation rates and environment 
JBalogh et all |2004| : iKauffmann ct al.; gOOJ) and via mea- 
surem ents out to z « 1 (jCucciati et al.l 120061 : ICooper et al.l 
I2OO7I ). 

Studies that investigated the relative importance of en- 
vironment on different scales indicated that small-scale en- 
vironment (Mpc scale) correlates str onger with galaxy prop- 
erties than large-scale environment (IKauffmann et al.l 2004 • 
Blanton et~al]i2006l: [sianton fc Berlindll2007l: IWilman etlD 



201(1 ). These results suggest that the mass of the dark 



matter halo in which a galaxy resides plays the dominant 
role in the galaxy's forma tion and evolution. For example, 
iBlanton fc Berlindl ((2003) found that the blue fraction of 
galaxies in rich groups and clusters is independent of envi- 
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ronment at fixed halo mass. In this paper we extend these 
results by focusing on galaxies in low-mass halos. 

The definition of environment is thus key to any dis- 
cussion of this topic. In this paper as in many other works, 
'environment' is defined as the local galaxy density at some 
smoothing scale. However, the physical interpretation of this 
definition of environment is not entirely straightforward be- 
cause galaxy density is a biased measure of dark matter 
density, in addition to being subject to shot noise and the 
effect s of redshift-space distortions (see, e.g., ICooper et al.l 
I2OO5I ). The choice of scale is also important. In this paper 
we make the distinction between (large-scale) environment 
and halo mass. We define environment as the local galaxy 
density on scales larger than the host dark matter halo of a 
given galaxy, using a fiducial smoothing scale of 10 h^^ Mpc. 
Halo mass is similar to 'small-scale environment', but the 
distinction is important. The dispersion of host halo masses 
for galaxies at a fixed 1 Mpc local density can be 1 dex or 
more and using distance to the n-th nearest neighbor yields 
an even larger spread in halo mass (jHaas et al.ll201in . This 
scatter demonstrates the need for quantifying halo mass ex- 
plicitly. 

The idea that galaxy properties depend primarily on 
host halo properties is physically motivated, since the 
halo virial radius corresponds to a physical transition be- 
tween infall and virialized motions, capable of support- 
ing strong shock front s and hot, thermalized gas (e.g., 
iDekel fc Birnboimll2006r ). In this sense, a galaxy group nec- 
essarily corresponds to a single, virialized dark matter halo. 
This is distinct from many other definitions of 'group envi- 
ronment'. For example, the 'Local Group' is not a group at 
all in this context, but rather a close pair of separate ha- 
los. The difference is an important one: for a halo to enclose 
both the Milky Way and Andromeda (with one galaxy at 
the halo center), it would have to be larger than lO'^'^ Moto 
equal their separation of ~ 800 kpc, a mass at which dif- 
fuse halo gas is heated to high temperatures, causing shock 
heating and stripping of the disks in our galaxy or M31 (or 
both). As a close pair of ~ 10^^ M© halos, their tidal in- 
teractions and gas physics is markedly different: there is no 
ram pressure and tidal stripping is much weaker. It is thus 
more informative to quantify halo mass rather than galaxy 
density on a small scale. 

However, environment does play a role in the formation 
of dark matter halos. In the numerical Universe, dark mat- 
ter halos exhibit 'assembly bias': at fixed mass, halo prop- 



erties depend on their formation histories (IGao et al. 



OL2 



Harker et al]|2006: Gao & White 2006; Wcchsler et all 



Wetzel et al.ll2007l : lLi et al.ll2008l : iDalal et al.ll2008l ). At low 



2005 



2006 : 



halo mass, older and more concentrated halos form in high 
density environments. At high mass the effect reverses, with 
younger, less concentrated halos forming in high-density re- 
gions, though the effect is also weaker. Semi-analytic models 
of galaxy formation predict that this assembly bias propa- 
gates into the evolution of galaxies at fixed halo mass, espe- 
cially in the formation histories and star formati on proper- 
ties of low-lum inosity galaxies in low-mass halos (|Zhu et al.l 
120061 : ICroton e t al. 200JZt)- ^^ first glance, these predictions 
appear in concert with the above observational evidence that 
the fraction of galaxies (at fixed luminosity) that are red and 
quenched is highest at high densities. 

Another approach to understanding the connection be- 



tween galaxies and dark ma tter is the Ha l o Occupation Dis- 
tribu t ion (HOD; see, e.g., Seliak 200d: Peacock fc SmithI 
200d: IScoc cimarro et al] I2OOII : iBerlind fc Weinberd |2002| : 
Coorav fc Shcth 20021 ). The HOD has emerged as our most 



powerful tool for understanding the clustering of galaxies 
and quantifying the bias between galaxies and dark matter, 
where bias here is defined as the ratio between the galaxy 
and matter densities on an arbitrary smoothing scale. The 
central quantity in the HOD is P{N\M), the probability 
that a halo of mass M contains TV galaxies within a de- 
fined sample. The galaxy sample may be arbitrarily defined, 
thus P{N\M) will be different for every sample definition. 
The bias of galaxies depends strongly on galaxy properties, 
with strong clustering for gala xies that are brighte r , red- 



der, and more ell i ptical (e.g.. J Norberg. et al.l I2OOII. 
Zehavi et all l2002l : IZehavi. et all l2005l : IZehavi et al.l 



2002 : 



201a : 



Li et al.ll2006l ). In the HOD context, these trends are all ex- 



plained by the different halos that different types of galaxies 
occupy. Brighter central galaxies live in more massive ha- 
los that are more highly clustered. At fixed galaxy mass, 
redder or more elliptical galaxies are more likely to be satel- 
lites in groups or clusters, thus enhancing their clustering. 
These observations fit into the HOD formalism without any 
explicit correlation between galaxy formation and environ- 
ment, other t han the fact that massive halos form at high 
densities (e.g. iBond et al.lll99ll ). 

The P{N\M) ansatz has been extended even fur- 
ther in the 'subhalo abundance matching' approach, 
which (in its simplest form) assumes a monotonic 
relationship b etween halo mass, M, and gala xy lu- 



minosity, L (iKravtsov et all l2004l: IWang et all |2006|: 



Vale fc Ostrikej I2OO6I: IConrov etlol I2OO6I: iMoster et all 



20ld: IConrov fc Wechslerl I2OO9I : IWetzel fc Whitd I2OI0I : 
Behroozi et al.ll20ld )rl Here, if the number density of ha- 
los above M is the same as that of galaxies above L, then 
M halos host galaxies of L brightness. They key to this 
approach is to include not just isolated halos, but also 'sub- 
halos', defined as bound, virialized halos located within the 
virial radii of larger host halos. The primary assumption of 
this approach is that (sub)halo mass entirely determines the 
luminosity of the galaxy within it, regardless of any other 
properties of the (sub)halo or the galaxy. 

The quantity P[N\M) makes the explicit, and strong, 
assumption that A'^ depends only on M. That is, the prop- 
erties of galaxies, including luminosity, color, age, and mor- 
phology, depend only on the mass of their host dark matter 
halo and have no additional dependence on a halo's for- 
mation history or environment. This formulation was not 
proposed with any theoretical prejudice about galaxy for- 
mation, but rather it was the simplest implementation and 
was successful in fitting the clustering properties of nearly all 
available data. However, the discovery of halo assembly bias, 
combined with the complex relation between galaxy prop- 
erties and environment, especially those governing galaxy 
growth and star formation , calls into question this ansatz. 

In [Tinker et al.l (|2008l '). we sought to test this P[N\M) 
ansatz through different galaxy clustering measures that 
probe different environments. We used the two-point cor- 



^ One can use halo maximum circular velocity rather than halo 
mass, or use galaxy stellar mass rather than galaxy luminosity. 
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relation function, ^g{r), to constrain halo occupation and 
make predictions for the void probability function. The pre- 
dictions were in excellent agreement with the measurements, 
demonstr ating the validit y of th e HO D from high to low 
densities. lAbbas fc ShethI (|2006l ') and ISkibba et all (120061 ) 
used other galaxy clustering measures to demonstrate that 
the assumption of P{N\M) is consistent with ob served data 
for lu mino sity-based ga l axy s amples. However, lYang et al.l 
l|2006l) and IWang et al] (|2008l ') found that the clustering of 
low-mass galaxy groups at fixed mass depends on the color 
of the galaxies contained within them such that groups with 
redder galaxies are more highly clustered. 

In this paper, we examine the relation between galaxy 
star formation histories, their host dark matter halos, and 
their large-scale environment. We explore in particular the 
extent to which the formation histories of galaxies relate to 
those of their host halos and whether the assembly bias seen 
in simulation extends to the galaxy population. We focus on 
low-mass galaxies, the regime in which halo assembly bias 
effects are the strongest, and we examine star formation his- 
tory based on narrow 4000-A break, D„4000, which, unlike 
color, is insensitive to dust reddening. In order to observa- 
tionally determine halo masses for SDSS galaxies and sepa- 
rate this from large-scal e environment, we construct galaxy 
group catalogs using the lYang et al.l (|2005r ) group-finding al- 
gorithm. Our use of a group finder is critical because galaxies 
of the same properties may occupy halos of a wide range of 
mass. At fixed stellar mass, galaxies may live at the center 
of low-mass halos or they may exist as satellites in group 
or cluster-sized halos. Anywhere between 10-40% of galax- 
ies a re satellites, depending on galaxy ma ss and color (see, 
e.g., lYang et al.ll2008l : IZehavi et all I2OI0I ) . The distinction 
between central and satellite galaxies will be the cornerstone 
of understanding the environmental correlations of galaxy 
properties. 

The evolution of satellite galaxies is more complex, 
and their properties also depend o n position within a 
halo (e.g., van den Bosch et al.l 120081 : [Hansen et al.l l2009l : 
IWeinmann et al.ll2010f ). Here we focus on the properties of 
central galaxies in halos below the group scale. We will ex- 
plore the observed dependence of satellite star formation 
rate on galaxy mass, halo mass, and halo-centric radius 
in detail in a companion paper (Wetzel, Tinker, & Con- 
roy 2011a), hereafter referred to as Paper II. Finally, in a 
third paper (Wetzel, Tinker, & Conroy 2011b), hereafter 
referred to as Paper HI, we will use these measurements 
combined with a high-resolution simulation to test various 
physical mechanisms for the quenching of satellite galaxies. 

We examine galaxies selected on either r-band magni- 
tude or stellar mass, M*, and we examine star formation 
history based on D„4000. Magnitude and D„4000 have the 
advantage of being well-defined observationally, though in 
Paper II we will focus entirely on stellar mass and instanta- 
neous specific star formation rate (SSFR), quantities which 
are derived and thus more model-dependent, but also more 
readily physically interpretable. We have examined all halo 
and environmental trends in both papers selecting on both 
Mr and M, as well as examining both D„4000 and SSFR. 
While these lead to slight quantitative differences in our re- 
sults, they do not change any of our results qualitatively. 

Throughout, we define a galaxy group as any set of 
galaxies that occupy a common dark matter halo (includ- 



Table 1. Volume-Limited Samples 
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Figure 2. The distribution of g — r color and D„4000 for Mr = 
[—19, —20] galaxies. The shaded region shows the distribution of 
color, with scale on the bottom x-axis. The red dashed curve 
shows the distribution of D„4000 for the same galaxies, with scale 
on the top x-axis. Even though 53% of galaxies are red {g — r > 
0.7), only 35% of galaxies are truly quenched (Dn4000 > 1.6). 
The blue and black dotted curves show the color distribution 
for quenched and active galaxies, respectively, according to their 
Dn4000 values. Quenched galaxies nearly always have red colors, 
while 23% of active g alaxies have red co lors, primarily because of 
dust extinction fe.e.. lMaller et al.llJOOsI 'l. 



ing cluster- mass halos), and we define a halo as having 
a mean interior density 200 times the background matter 
density. A host halo is a halo that is distinct: its center 
does not reside within the radius of a larger halo. A sub- 
halo is one whose center is located within the radius of a 
larger halo. For all calculations we assume a flat, ACDM cos- 
mology of (n„, (78, ilb, Us, ho) = (0.27, 0.82, 0.045, 0.95, 0.7). 
For galaxy magnitudes we do not assume a value of h, 
but for brevity we write all magnitudes as A'U rather than 
Air — 51og/io. Stellar masses are in units of Mq/H'^. 



2 DATA AND MEASUREMENTS 

2.1 NYU Value-Added Galaxy Catalog 

To construct our galaxy sample s, we use the NYU Value- 
Added Galaxy Catalog (VAGC: iBlanton et al] 120051 ) based 
on the spectroscopic sample in Data Release 7 (DR7) 
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Figure 1. Left panel: R elation between r-band magnitude, Mr, and stellar mass, M, , with the latter taken from the kcorrect code of 
iBlanton fc Roweij 1120071 ). Vertical lines show the magnitude limits for our volume-limited samples as listed in Table 1, while horizontal 
lines indicate the stellar mass limits within each volume-limited sample. Right panel: Cumulative numb er density of galaxies a s a function 
of both Mr and M*, offering a rough conversion between the two. The star indicates M* from the iBlanton. et al.l 1 12003) luminosity 
function. 
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Figure 3. Distribution of Dn4000 values for galaxies in mag- 
nitude bins. Faint galaxies arc predominantly active while bright 
galaxies are predominantly quenched but still have a significant 
active population. All sample distributions are bimodal with a 
minimum at Dn4000 ~ 1.6. We refer to galaxies with D„4000 > 
1.6 as 'quenched'. 



of th e Sloan Digital Sky Survey (SDSS; lAbazaiian et al.l 
I2OO9I ). Specifically we use dr72bight34. We construct three 
volume-limited samples, as listed in Table 1, which con- 



Figure 4. Quenched fraction, /q, defined by Dn4000 > 1.6, 
as a function of large-scale (10 h~^ Mpc) density of galaxies. At 
p/p > 1, /q rises monotonically with density, while in underdense 
regions the quenched fraction is nearly constant. 



tain all galaxies brighter than Air = —18, Mr = — 19 and 
Mr ~ —20, respectively. Within each volume-limited sam- 
ple, we construct another sample that is complete in stellar 
mass. The stellar masses are also taken fr om the VAGC and 
are derived from the kcorrect code of IBlanton fc Roweid 
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((20071), which assumes a IChabrieiJ (|2003l ) initial mass func- 
tion. Fig. [l] (left) shows the comparison between stellar mass 
and r-band magnitude. The vertical and horizontal lines in- 
dicate the magnitude and stellar mass limits, respectively, 
for our samples. By constructing our stellar mass samples 
this way, they contain significantly fewer galaxies than the 
luminosity-defined samples. This limits the statistics and ef- 
fective halo mass range of the stellar mass samples, but this 
method has the benefit of requiring no weighting scheme 
for red galaxies that are fainter than blue galaxies at fixed 
stellar mass. We will construct groups and present results 
based upon both sample definitions in order to demonstrate 
robustness to galaxy property selection. To allow convenient, 
though approximate, conversion between the two properties. 
Fig. [T] (right) also shows the cumulative number densities of 
galaxies based upon Mr and M*. 

For galaxy pairs that are too close to obtain spectra 
because of the 55 arcsecond width of SDSS fibers ('fiber 
collisions'), we use the internal correction to the fiber cor- 
rections within the VAGC, namely that the collided object 
is given the redshift of the nearest galaxy in terms of an- 
gular separation, provided that this redshift is in agreement 
with the ph otometric redshift ob tained by with the SDSS 
photometry (jBlanton et al.ll2005l ). 

2.2 Bimodality of Galaxy Properties 

Galaxies can be roughly divided into two distinct color 
categories: the 'blue cloud' and the 'red sequence'. The 
former is comprised of galaxies that are actively form- 
ing stars at the current epoch, giving them blue colors. 
These galaxies are primarily disk-dominated with signifi- 
cant amounts of cool gas. The red sequence is comprised 
(mostly) of galaxies with old stellar populations, devoid 
of cold gas, and usually exhibiting elliptical morphology. 
The bimodality of the galaxy color distribution is well mea- 
sured at z = jS tratcva, et al.ll200ll: 'Blanto n. et al.ll2003l : 
iKauffmann et al.l l2003: M adgwick et a l. 2003) and has been 
shown to exist at z = 1 (JBeU et al.ll2004l : ICooper et al.ll2006l: 
Willmer. et al .1120061 ') and even up to z = 2 (jWilliams et al.l 



2009[ ). 



Using galaxy color as a proxy for star formation ac- 
tivity can be problematic, as dust reddening can cause a 
:as-rich disk galaxy to be classified as a red sequence object 
Mailer et al.ll20od : [Masters et al.ll2010l ). To avoid this prob- 



lem, we use D„4000, which is derived from SDSS spectra. 
D„4000 is a diagnostic of the light- weighted age of the stellar 
population and thus is sensitive to the integrated star forma- 
tion history of the galaxy. We obt ain this quantity from th e 
JHU-MPA spectral reduction^ (|Brinchmann et al.1 IJOol ). 
Fig. [2] compares the distribution oi g — r color to that of 
D„4000 for —19 > Mr > —20 galaxies, our fiducial magni- 
tude bin throughout this paper. The red sequence is defined 
a strong peak in the color distribution at g — r ~ 0.9 and at 
D„4000 « 1.9. Both distributions are clearly bimodal, but 
based on color alone one would conclude that the majority 
of galaxies in this magnitude range are quenched. The pic- 
ture painted by D„4000 is strikingly difi'erent: the majority 
of galaxies of this magnitude are instead active, but many 



^ |http : //mm . mpa-garchlng ■ mpg . de/SDSS/DR7/ 1 
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have been dust-reddened. We have confirmed that this popu- 
lation of active but red galaxies is caused primarily by dusty 
spirals by applying Galaxy Zoo morphological classifications 
(JLintott et al.ll201ll ) to the sample: ~ 70% of are identified 
as s pirals, and ~ 50% a re identified as edge-on spirals (see 
also [Masters et al.l bOlOD . 

One concern in using galaxy properties defined by spec- 
tra is fiber collisions, which cause ~ 8% of all galaxies in the 
survey to not have spectra. To correct for this, for each fiber- 
collided galaxy we draw a random value of D„4000 based on 
a sample of spectroscopic galaxies that are within ±0.25 in 
Mr and ±0.05 in p — r of the collided object. This method 
allows us to use the colors to infer statistically the spectro- 
scopic properties of the collided sample. This method refiects 
the fact that not all collided red galax;ies are truly quenched, 
and it preserves the fact that the color distribution of col- 
lided objects (and the distribution of D„4000) is distinct 
from the overall color distribution at the same magnitude. 
Using this method, our results are generally consistent with 
those obtained by simply removing the collided galaxies from 
the samples after constructing the group catalog. 

Fig. [3] shows D„4000 distribution for galaxies of differ- 
ent absolute magnitudes. As seen in previous studies, the 
bimodal distribution in sub-L* galaxies is apparent but is 
dominated by an younger, active population. At L > L,, 
the bimodality is still present but has been tilted to the 
older, quenched population. The brightest galaxies are dom- 
inated by an red, old population that is not forming stars 
at any appreciable rate. Because the minimum of the bi- 
modality occurs at D„4000 = 1.6 across all magnitudes we 
explore, we define galaxies with D„4000 < 1.6 as 'active' 
and Dn4000 > 1.6 as 'quenched'. 

Another possible concern with galaxy properties defined 
by SDSS spectra is the finite 3 arcsec aperture of the fibers, 
corresponding to ~ 2/i~^kpc at the redshifts we consider. 
Thus, for many galaxies the SDSS fiber covers only the inner 
light profile, where the stellar population may by older (or 
younger) than the rest of the galaxy. As can be seen in Fig.|31 
the fraction of blue galaxies that are labeled as quenched via 
D„4000 is negligible compared to the fraction of red galaxies 
that have active stellar populations. Thus, while aperture 
effects may bias D„4000 values somewhat, the effect of dust 
obscuration on color is far more important. 



2.3 Measuring Large-scale Environment 

For each galaxy, we determine the large-scale environment 
by counting the number of neighboring galaxies within a 
sphere of radius 10 h~^ Mpc centered on each galaxy. This 
quantity is a biased indicator of the dark matter density 
field, but at 10 h~^ Mpc this bias is a simple linear factor 
and any stochasticity is minimal. We count the number of 
galaxies above the corresponding magnitude threshold for 
the each sample, and so the tracer of the density field has a 
different bias for each sample. We do not correct for this be- 
tween the samples, but note that the relative bias between 
Mr < —18, —19, an d —20 samples is at the ~ 5% level 
(jSwanson et al.ll2008l ). This galaxy density measurement is 
affected by galaxy peculiar velocities, but this effect is min- 
imal at 10 h~^ Mpc, as we demonstrate in Appendix A. 

We also choose 10 h~^ Mpc as it represents a clear dis- 
tinction from a galaxy's small-scale environment as encapsu- 
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Figure 5. Quenched fraction of satellite galaxies, /rf*, as a func- 
tion of A/host for bins in galaxy magnitude. At fixed halo mass, 
brighter satellites are more likely to be quenched, and satellites of 
all magnitudes show a similar increase in quenched fraction with 
halo mass. 



lated by its host halo, especially since 10 h~^ Mpc represents 
the maximum extent from which the most massive galaxy 
clusters have accreted their mass. In tests we find that our 
results show little dependence on the exact smoothing scale 
chosen. In Paper II we will further discuss the properties of 
galaxies on the outskirts of groups and clusters. 

To correct for survey geometry and incompleteness, we 
use of random catalogs. For each volume-limited sample, we 
produce a catalog of lO'^ random points distributed with the 
angular selection function of SDSS DR7 using the angular 
mask provided with th e VAGC in combinatio n with the soft- 
ware package mangle (jSwanson et al.ll2008l V Each random 
point is also assigned a random redshift such that the co- 
moving space density of randoms is constant with redshift. 
For each galaxy, we correct for incompleteness by multi- 
plying the observed number of galaxies by the ratio of the 
number of random points divided by the expected number 
of randoms if the completeness were unity. The large num- 
ber of random points ensures that shot noise within each 10 
h^^ Mpc sphere is at the sub-percent level. 

Fig-Hlshows the fraction of quenched galaxies as a func- 
tion of large-scale density. For all magnitude bins there is 
clear dependence on density such that quenched galaxies 
preferentially live in overdense environments. At p/p > 1, 
the slope of the quenched fraction with log p/p is roughly 
independent of galaxy magnitude. At p/p > 10, the frac- 
tion of quenched galaxies in the faintest magnitude bin rises 
steeply. Other studies have also noted that the faintest red 
galajci es predominantly are found in the hig hest density re- 
gions (|Hogg et aI.ll2004l : lBlanton et aI.ll2005lV 



3 GROUP FINDING ALGORITHM 

Motivated by the framework of the HOD model, we assign 
ha los to galax i es usin g the group-finding algorithm detailed 
in lYang et al. (120051) a nd applied to the SDSS in a series 
of papers jYang et all 1 2 007. 2008. 2009: Weinmann etld] 
I2OO6I : IWang et all bOoSTvan den Bosch et ~l.ll2008l ). Giwn 
the known statistics of dark matter halos, the halo masses 
of groups of galaxies can be retrofit to the observed distribu- 
tion of galaxies. Because this algorithm has been described 
in detail in the papers above, in this section we give a syn- 
opsis of the group finder and the minor modifications we 
made to the algorithm, and we provide a more thorough 
description in Appendix B. 

Our s t arting point differs from the algorithm of 
lYang et al.l (|2005l ). To begin, we ascribe to each galaxy in 
the VAGC a dark matter halo mass via the subhalo abun- 
dance matching method described in §1. Typically, abun- 
dance matching is used to put galaxies into a population of 
dark matter (sub)halos, but we invert this process to put 
(sub)halos around galaxies as the first step in the group- 
finding algorithm. The abundance matching method relates 
galaxy luminosity (or stellar mass) to (sub)halo mass by 
assuming that the two are monotonically related. Thus, a 
(sub)halo of mass Mo contains a galaxy of luminosity Lq 
such that 



^(L) dL - 



ntot{M)dM 



(1) 



M„ 



where $(!/) is the galaxy luminosity function (or alternately 
the stellar mass function) and ntot(M) is the halo plus sub- 
halo mass function. 

In applying equation ([1} we use the lTinker et al.l (|2008l ) 
function fo r the halo mass functio n and the subhalo mass 
function in [Tinker fc Wetzel (|2010l ) (their equation 12; but 
we have changed the normalization from 0.3 to 0.2 to ac- 
count for the differences in the subhalo mass function be- 
tween redshift 1 and 0). The total number density of halos 
of mass M is then 

ntot(M) = n(Mhost) + / n(Afsub|Mhost)n(Afhost) dMhost (2) 

where we have denoted halos that are distinct as 'host' ha- 
los and halos that are contained within host halos as 'sub'. 
Thus, each galaxy is given a halo mass regardless of whether 
it is a satellite within a larger host halo or resides in the 
fielcEI. 

To obtain the left-hand side of equation ((T} , we sum all 
galaxies, Li, brighter than Lq, weighed by Knax to obtain 
the number density, given by. 



^{>Lo)^ J2 Vn.UU) 



(3) 



Li>=Lo 



where Knax is the maximum volume to which each galaxy 
is observed, as calculated in the VAGC to account for sur- 
vey incompleteness. Once again, luminosity can be replaced 
by stellar mass, and a comparison of the cumulative den- 
sity of galaxies as a function of magnitude and M» is given 



^ Note that M^uj, is the mass of the subhal o at the time of ac- 
cretion. See details in lTinker &: Wetzell 1I2OIOI ) and a discussion of 
subhalo finding in Paper III. 
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Figure 6. Quenched fraction, /q, vs. 10 h~^ Mpc ovcrdcnsity, p/p, for various bins in galaxy magnitude. In each panel, the solid curve 
is the overall /q-p/p relation, with the filled circles showing /q°" and /q^'. Bright central galaxies do show a slight increase in quenched 
fraction with density. 



in Fig. [T] Although we later cull the VAGC to create the 
volume-liniited samples in Table 1, using the full flux-limited 
catalog suppresses sample variance in the calculation of the 
luminosity function for the initial halo mass assignment. 

In the lYang et al.l (|2005J) methodology, this initial 
matching of galaxies to groups and halo mass is done with 
a redshi ft-space friends-of-fr iends (FOF) linking algorithm 
(see also iBerlind. et al.ll2006l V Using the inverse abundance 
matching approach, we flnd results that are consistent with 
the FOF algorithm. See the tests in Appendix C. 

Once a halo has been assigned to each galaxy, each 
galaxy has an associated halo mass, virial radius, and ve- 
locity dispersion via the virial theorem. We then determine 
the probability that each galaxy is a central galaxy in a 
host halo or a satellite galaxy in a subhalo. If, projected 



on the sky, a galaxy lies within a more massive galaxy's 
virial radius, we determine an angular probability that the 
galaxy is a satellite by assuming that the number density 
profile of sat ellite galaxies follows the dark matter given 
by an NFW (JNavarro et al.lll997bh pr ofile. We assume the 
concentration-mass relation given by iMaccio et al.l pOOST l 
for our cosmology, but note that the results are insensitive 
to this choice. We also assign a line-of-sight satellite prob- 
ability to the galaxy given its redshift offset from the more 
massive galaxy, where we assume that the host halo's satel- 
lites are distributed in a Gaussian along the line of sight. 
If the product of the angular and line-of-sight probability is 
above a calibrated constant, then the galaxy is considered to 
be a satellite in the larger host halo. Once we have applied 
this routine to all galaxies and we have a list of candidate 
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Figure 7. Left panel: The relation between quenched fraction, /q, and 10 h^^Mpc ovordensity, p/p, for Mr = [—19,-20] galaxies. 
The solid curve shows the overall relation from Fig. |4] The different filled circles show the quenched fraction divided into central and 
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(this correction is described in detail in Appendix C). Right panel: Same as left panel, but using the galaxy group catalog defined with 
stellar mass as opposed to galaxy luminosity as the tracer of group mass. The number density of galaxies with A4t > IO^^Mq//)^ 



~ 30% lower than that of galaxies with M, 
mass halos. 



< — 19 (see Fig.[T]|, thus the central galaxies in the right hand panel occupy slightly higher 



host halos, we recompute host halo masses, again by the 
abundance matching method. However, from this point on, 
we assume instead that host halo mass correlates with the 
total luminosity of the group. Thus, we use the halo mass 
function in equation lU for host halos only (no subhalos), 
and the luminosity now is the total luminosity of the galax- 
ies within the halo. We iterate this procedure until host halo 
masses for the groups have converged. Our catalog contains 
groups varying in mass from 2 x 10^^ to 10^^ H^^Mq. 

For each sample, we perform the group-finding proce- 
dure once using total galaxy luminosity as the tracer of halo 
mass and once using stellar mass as the tracer of halo mass. 
There is evidence that stellar mass correlates stronger with 
halo mass than l uminosity, especially when dividing a galax y 
sample by color (|Mandelbaum et alJl2006l : lMore et al.ll2010l ). 
However, there are also significant uncertainties in determin- 
ing stellar mass fr o m broadband optical magnitudes (e.g., 
IConrov et ai]|2009l : IConrov fc GunnlbOlOl '). and so it is not 
clear whether stellar mass or luminosity is preferable here. 
Nonetheless, we will demonstrate that all our conclusions 
are robust to choice of halo mass tracer. 

Note that our algorithm defines the brightest (or most 
massive) group galaxy as the central galaxy, that is, the 
galaxy residing at the minimum of the host halo's potential 
well. For halos with Mhost < 10^^ /i~^Mq, central galax- 
ies are expected to be the brightest galaxy in the halo the 
vast majority of th e time (^Zchavi, ct al. 2 0051 : [Tinker et al.l 
l2005l : IZheng et al.i r2005: Yang ct al. 2009). For cluster-sized 
halos, Mhost ^ 10^* h^^ M©, wheth er the brightest ga laxy 
is always central galaxy is not clear (jSkibba et al.ll201ll ). In 
statistical terms, halos of this mass are rare and do not con- 
tribute significantly to the results presented here. Because 



we focus on central galaxies with L < L* luminosities, this 
is not a concern for our results. 



lYang et al.l (|2005l ') and subsequent papers have demon- 
strated the efficacy of the halo-based group finder, but the 
approach is not infallible. The algorithm is designed to ob- 
tain the correct mean number of galaxies per halo, but be- 
cause of projection effects and redshift-space distortions, 
some galaxies will erroneously be labeled as 'central' when 
in fact they are within the virial radius of a larger halo, 
while other central galaxies will falsely be labeled as 'satel- 
lites'. This misclassification of galaxy type means that the 
measured fraction of central and satellite galaxies that are 
quenched will be biased because satellite galaxies are more 
likely to be quenched than central galaxies at fixed M*. In 
Appendix C we use a high-resolution A*'-body simulation to 
quantify the false-classification rate and demonstrate how to 
statistically correct the measurements of /q for this effect. 
On average, ~ 10% of galaxies are misclassified, and all of 
the results in this paper include a correction for this misclas- 
sification, except where stated otherwise. We note that the 
group-finding algorithm is actually more efficient at iden- 
tifying the central galaxies than satellite galaxies within a 
given sample, yielding a completeness of ~ 95% for most 
central galaxies (see App. C). 

We have also checked our r esults against a SDSS DR7 
version of the I Yang et al.l (|2007l ) group catalog (kindly sup- 
plied by Frank van den Bosch), incorporating the same 
kcorrect stellar masses D„4000 measurements we use here 
but retaining differences in assumed cosmology, halo virial 
definitions, and group-making methodology. All of the re- 
sults are consistent within errors, demonstrating the robust- 
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Figure 8. Quenched fraction of satellites, /q**, as a function of 
10 h^^ Mpc ovcrdensity, p/p, for satellites with Mr = [—19, —20]. 
Filled circles show the measurements from Fig. [T] Solid curves 
show /q^ in bins of log M^ost ■ 



ness of the iterative group-making method to the details of 
the initialization of the procedure. 



4 RESULTS 

4.1 Properties of Satellite Galaxies at Fixed Halo 
Mass 

For our fiducial results we examine a bin in magnitude of 
Mr — [—19,-20]. These galaxies are faint enough that they 
span a wide range of host halo masses (as satellite galaxies) 
but bright enough that a volume-limited sample contains 
sufficient statistics for fine binning in halo mass. By looking 
at galaxies in a relatively narrow range of luminosities, we 
are restricting our analysis to galaxies that are in a narrow 
range of M^ub (for satellite galaxies) and Afhost (for central 
galaxies) . 

Fig. [5] shows the fraction of quenched satellites, /q'', 
as a function of A/host for bins in AU- The slope of /q"' 
is independent of galaxy magnitude, but the amplitude in- 
creases monotonically with Mr. These results indicate that 
accretion of a galaxy onto a larger halo contributes signifi- 
cantly to the buildup of quenche d, passive galaxies, in agree- 
ment with pre vious works (e.g.. Ivan den Bosch et al.l 120081 : 
[Tinker fc Wetz el ,201Q) . 

The satellite trends in Fig. [S] are broadly consistent 
with previous results on galaxy mass and halo mass de- 
pende nce. Using the maxBCG cluster sample. [Hansen et al.l 
l|2009l) found that the red fraction of galaxies increases 
with cluster richness, though their red fraction is somewhat 
higher than our quenched fraction due to the use of color 
rather than a dust-insensitive diagnostic. Other works using 
more direct star formation rate indicat o rs have found more 
similar trends jWeinmann et all l2006l : iKimm et al.l |2009| : 



Figure 9. Solid and dashed curves indicate the mean logM^ost 
as a function of 10 h~^ Mpc overdensity, p/p, for satellite galax- 
ies and all galaxies in the Mr = [-19,-20] bin, respectively. 
Here, (logMi^ost) is a galaxy- number weighted mean. Unlike cen- 
tral galaxies, the typical halo mass probed by satellite galaxies 
increases with density. The filled circles indicate the fraction of 
all galaxies that are central galaxies as a function of p/p, as given 
by the right-hand y-axis. 



IWeinmann et al.ll2010l : Ivon der Linden et al.ll2010l '). We will 
examine trends of satellite quenching in much more detail in 
Papers 11 and IIL For now, these results will be important 
when interpreting the correlations with density in Fig. |4l 



4.2 Dissecting the Correlations with Density 

Fig. |5] shows the /q-p/p relation, broken into central and 
satellite galaxies, for four magnitude bins. For L < L* galax- 
ies, the quenched fraction of central galaxies is independent 
of large-scale density, spanning the entire range of environ- 
ments from the deepest voids (p/p ~ 0.1) to cluster infall 
regions (p/p > 10). The entire correlation with environment 
is driven by the satellite galaxies. The results are consistent 
with the scenario in which all environmental correlations are 
due to the change in the halo mass function. However, for 
the Mr = [-20.5,-21] magnitude bin there is a clear de- 
pendence of /q°" with density, increasing from /q = 0.55 to 
0.65 over two decades in density. We will discuss the brighter 
galaxies subsequently, focusing now on understanding the 
/q ~ p/p correlation for our fiducial sample. 

The left panel in Fig. [7] shows the /q-p/p relation for 
Mr = [-19,-20] galaxies taken from Fig. [Hb. The dotted 
curves indicate the raw quenched fractions, uncorrected for 
satellite-central mislabeling. Because mislabeling only in- 
creases /q°" and decreases /q'', these curves can be con- 
sidered upper and lower limits on these two quantities, re- 
spectively. The right panel shows the same breakdown for 
groups defined by stellar mass rather than luminosity. In 
both panels, the dependence of /q on environment is caused 
entirely satellite galaxies; /q** rises rapidly when p/p > 1. 
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Understanding this satellite dependence on p/p requires 
a decomposition into host halo mass. As was demonstrated 
in Fig. O the quenched fraction of satellite galaxies at fixed 
luminosity is a strong function of Mhost. Fig. [8] breaks /q'* 
into bins of fixed Mhost. When controlling for halo mass, the 
quenched fraction of satellites also is independent of envi- 
ronment within the statistics of the sample. Thus, the ob- 
served correlation between /q and p/p simply arises from 
the change in the halo mass function with environment. 

Fig. [9] shows the relative number of centrals and satel- 
lites as a function of large-scale environment. Unsurpris- 
ingly, central galaxies dominate the statistics in the voids. As 
p/p increases, the central fraction monotonically decreases, 
dropping under 50% at p/p > 20. Because we are look- 
ing at a bin in galaxy magnitude, the halo mass probed by 
centrals is constant with density at logAf ~ 11.7, but the 
typical halo mass probed by satellite galaxies monotonically 
increases with density. At 10 times the mean galaxy density, 
(log Mhost) ~ 13 for satellite galaxies. 

The change in (log Mhost) with large-scale density is 
a well-understood phenomenon: more massive dark mat- 
ter halos are mor e likely form in higher o verdensities (e.g., 
iBond et all Il99ll : ISheth fc TormenI |2002| V Fig. [5] demon- 
strates that /q depends on Mhost, thus as p increases the 
quenched fraction of satellites increases as a result of the 
changing halo mass function. 
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Figure 10. Quenched fraction of bright central galaxies vs. 
10 h~^ Mpc overdensity. The circles represent f^'^ for all cen- 
tral galaxies in this magnitude bin, while the curves show these 
galaxies split into different halo masses, demonstrating that the 
increase in quenched fraction with density also persists for fixed 
halo mass. 



4.3 Bright Central Galaxies 

As demonstrated in Fig. [6] the quenched fraction of L > L, 
central galaxies does exhibit a weak but significant corre- 
lation with density. Using our Poisson error bars, a line 
with slope d/q'Vrflogp = 0.08 yields a Ax^ of 25 relative to 
the best-fit constant value. Due to scatter between Mr and 
Mhaio, central galaxies in this magnitude bin span a range 
in halo mass. Thus we must determine if this correlation is 
due to an intrinsic relationship between /q and density or a 
change in the typical halo mass probed in low and high den- 
sities. Fig. [TOl shows the results for the Mr = [-20.5,-21] 
central galaxies broken into bins of fixed halo mass. For each 
halo mass bin, the quenched fraction shows a correlation 
with density comparable to that of the overall magnitude 
bin, which demonstrates that more massive galaxies exhibit 
a true correlation with environment at fixed halo mass. 



5 CONFRONTING THEORY WITH 
OBSERVATIONS 

5.1 Quenched Galaxy Fractions vs. Old Halo 
Fractions 

Although a full model of galaxy formation requires incorpo- 
ration of many unsolved physical processes, it is possible to 
test the ansatz that the oldest (or slowest growing) galaxies 
reside in the oldest halos. In the context of the observational 
metric used here, 'old' galaxies have ceased growing in mass 
through star formatior|j, while 'old' halos have accumulated 



* Low-mass, quenched galaxies may increase their mass through 
merging, but signific ant mass growth through mergers is un likely 
for sub-L, galaxies l lMallei]l20o"8l : IConrov fc WechsleJlJOOgl) and 



most of their mass at early times and have little present 
mass growth. 

In this section we compare theoretical predictions to 
results from the group catalog that uses stellar mass rather 
than luminosity to estimate halo mass. We create predic- 
tions using halos in a high-resolution A'^-body simulation 
described in Appendix C and Paper III. This simulation has 
sufficient mass and spatial resolution to track the evolution 
of halos and subhalos down to 10^^ /i~^MQ. We populate 
the halos and subhalos of this simulation with galaxies using 
the abundance matching tec hnique described in jj2, assigning 
galaxy luminosities with the lBlanton. et al.l (|2003h luminos- 
ity function. This allows us to calculate p/p for objects in 
the simulation in the same manner as our DR7 resultqj. 

Fig. Illl shows fQ-p/p for three bins in galaxy mass (all 
less massive than the knee in the stellar mass function). In 
each stellar mass bin, the quenched fraction is independent 
of environment, consistent with the results from from §4. 
From the group catalog, these galaxies live in halos of mass 
log Mhost = 11.4, 11.5, and 11.6. We bin the halos in the 
simulation to match the halo masses probed in the stellar 
mass samples in Fig. 1111 In each bin, the halos are divided 
into 'old' and 'young' using the fractional growth rate of 
each halo. In each panel, we plot the old fraction of halos, 
/old, as a function of p/p, normalized such that the overall 
old fraction of halos in each mass bin is the same as /q"" of 



major mergers would more likely involve a gas-rich galaxy, leading 
to significant star formation. 

^ When measuring p for galaxies in the stellar mass limited group 
sample, we use the full volume-limited galaxy sample (to the mag- 
nitude limit of the volume) in order to maximize the statistics 
with which we measure density. 
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Figure 11. Central quenched fraction, /q"" vs. 10 h ^ Mpc overdensity, p/p, in three bins of stellar mass. Filled circles indicate the 
measurements from the group catalog. The mean logarithmic halo mass for each stellar mass bin is 11.4, 11.5 and 11.6 for panels (a), 
(b), and (c), respectively. Curves indicate the 'old' fraction of dark matter halos in the corresponding mass ranges, where 'old halos are 
defined as those with the lowest growth since a given redshift, and the growth threshold is chosen to match the corresponding fraction 
of quenched central galaxies. The solid curves show halos with lowest growth from z = 0.5 to z = 0, while the dotted and dashed curves 
are from z = 1.0 and z = 0.3, respectively. For halos, p/p is measured in the same manner as the data, using galaxy density on the 10 
h~^ Mpc scale. Galaxies arc placed into halos and subhalos using Eq. ((ij. All errors are Poisson. Galaxies with the lowest growth rates 
do not simply reside in halos with lowest growth rates. 



the galaxies. In each panel, the /om depends strongly on p/p 
such that older halos are preferentially found in high-density 
environments. This is consistent with earlier numerical re- 
sults on assembly bias in which old, low-mass halos are more 
highly clustered than young halos of the same mass (see the 
references in §1). The different curves in each panel represent 
different redshift baselines over which to calculate the halo 
growth rate, definecQ as 1 - M{z)/M{0), where z = 0.3, 
0.5, and 1.0. Using a shorter redshift baselines leads to a 
somewhat stronger environmental dependence. But regard- 
less of the definition of halo growth rate, the trend for the 
halos is not matched by the galaxies. This implies that halo 
formation history is not correlated with galaxy formation 
history; at fixed stellar mass, old, quenched galaxies and 
young, star-forming galaxies live in halos drawn effectively 
randomly from the distribution of halo growth rates. 

To demonstrate this result more clearly. Fig. [12] shows 
the ratio of the quenched fractions in high- and low-density 
regions as a function of central galaxy stellar mass. The 
ratio for SDSS centrals is near unity at all Af*. For galaxies 
logM» < 10, the results are inconsistent with halo growth 
being correlated with galaxy growth. 



^ Note that the halo mass used, M{z), is the maximum mass 
the halo has achieved up to redshift z. The maximal mass should 
correlate closer to gala xy stellar mass than instantaneous mass 
dWetzel fc Whitg|2010l) . Halos that lose mass from redshift z to 
zero would then have growth rates of 0. The rank-ordering of halos 
by growth rate is insensitive to choice of instantaneous versus 
maximum mass. 



5.2 Halo Formation Times Versus the Redshift of 
Galaxy Quenching 

Rather than assuming a monotonic relationship between 
galaxy age and halo age, we can also test whether halo for- 
mation times correlate with inferred galaxy quenching times 
for central galaxies. We define quenching time as the look- 
back time to the redshift at which a galaxy last formed stars. 
This comparison is model-dependent, both on the observa- 
tional and theoretical sides, but we will use quantities and 
parameters to maximize the validity of the comparison. 

To model quenched galaxies with stellar population syn- 
thesis (SPS) models, we assume a single burst of sta r for- 
mation. Using the SPS models of IConrov et al.l ((20091) and 
IConrov fc GunnI (|20ld ) we calculate the time evolution of 
D„4000. F or these c a lculat ions, we assume an initial mass 
function of lChabrieJ (|2003l ) and solar metallicity. To model 
actively star forming galaxies, we assume a star formation 
history of SFR ex t x ex p(— f/r) wher e r = 3.5 Gyr, consis- 
tent with the results of iNoeske et al.l (J2007) for the stellar 
masses probed here (see Paper III for full details regarding 
the determination of r). Once again, we assume solar metal- 
licity in order to calculate D„4000 as a function time. With 
these models, we create a lookup table to determine either 
the quenching timescale or formation timescale for passive 
and active galaxies, respectively. These timescales are model 
dependent, being especially sensitive to the metallicity as- 
sumed. Thus, the quantity we compare to theoretical pre- 
dictions is the variation relative to the mean timescale as 
a function of density. This quantity removes a significant 
portion of the model dependence. 

For dark matter halos, a common quantity used to mea- 
sure the formation time of a halo is ii/2 the time since a halo 
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accumulated half its current mass. However, for this mea- 
sure the peak in the age distribution is ~ 9.5 Gyr {z ~ 1.6). 
Most low-mass qu enched galaxies ha ve migrated to the red 
sequence aX z < 1 (|Drorv et al.ll2009l ') , demanding a measure 
of halo formation time that probes more recent structure 
formation. To this end, we use fss the time (in Gyr) since 
a halo reached 85% of its maximal mass. For this measure, 
the peak in the halo age distribution is ~ 4 Gyr {z ~ 0.4), 
but with a broad tail to higher redshift. 

Fig. 113b shows the mean D„4000 values as a function of 
p/p for low-mass quenched central galaxies (D„4000 > 1.6). 
The mean D„4000 does not vary with density, implying that 
the quenching time is constant with density as well. The 
mean quenching timescale is constant to within ~ 100 Myr 
across all densities. Fig. 113b compares these results to the 
formation times of dark matter halos. The solid curve is 
the change in fgs relative to the mean for all halos with 
logMhaio = 11-5 ± 0.15 (equivalent to the host halo masses 
for galaxies with M* = IO^'^Mq/Zi'^). In this comparison, 
the key assumption is not that the oldest halos contain 
quenched galaxies, but rather that the quenching time cor- 
relates with halo formation epoch. Under this assumption, 
the mean quenching timescale should vary by roughly 1.3 
Gyr with density. If we combine these two assumptions we 
obtain the dashed curve — here, we isolate the oldest 20% of 
galaxies in each bin in p and determine Atgs for only these 
halos. In this scenario, the mean formation time varies by 
nearly 3 Gyr from low to high densities. 

The dotted curve plots Afi/2 rather than using tgs, 
which shows no correlation with density. However, as we ar- 
gue above, using a timescale that probes structure at 2 > 1 
is unlikely to have any correlation with the recent star for- 
mation and growth of galaxies, in contrast to our chosen 
statistic of igs . Recent halo growth correlates strongly with 
current environment, but if one measures halo growth so far 
in the past, these correlations become washed out unless one 
probes the wings of the distribution. 

Fig. 113b is parallel to ll3t i. but now for star-forming cen- 
tral galaxies (D„4000 < 1.6). Once again, the mean D„4000 
is nearly independent of density. When converting to for- 
mation timescale, however, a slight positive trend with p 
is apparent; the timescale varies by roughly 0.6 Gyr from 
low to high densities. Results for dark matter halos show a 
stronger trend with density, especially at p/p > 10. Under 
the assumptions of our stellar population synthesis mod- 
els, star forming galaxies in high densities are slightly older 
than their counterparts in underdensities, but the amplitude 
of the effect is about half what one would predict using the 
formation times of dark matter halos. 



5.3 Evolution of the Stellar-to-Halo Mass Ratio 

Figs. [TT] and [13] imply that halo growth rates are not corre- 
lated with galaxy growth rates for galaxies below the knee 
in the stellar mass function — or, perhaps more straightfor- 
wardly, that the mechanisms the quench galaxy star for- 
mation are not correlated with dark matter growth. If halo 
growth rate is indeed uncorrelated with galaxy growth rate, 
then the M«/Mhaio ratio for low-mass galaxies will evolve 
differently for central quenched galaxies as compared with 
those that remain active. For quenched galaxies, M, is es- 
sentially fixed while Mhaio continues to evolve, while ac- 



tive galaxies continue to grow at a mu ch higher fractional 
rate t han their halos (see, e.g.. Fig. 4 in lConrov fc Wechslen 
l200g ). Assuming that a halo is on the mean M«/Mhaio re- 
lation at the time its star formation is quenched, Fig. 1141 
shows how far off the mean Mt/Mhaio relation at z = that 
halo would be as a function of the redshift of quenching. 
The details of this calculation are as follows: first, we set up 
the M,/Mhaio ratio as a function of halo mass through Eq. 
lU , using the halo plus subhalo mass func tion in Eq. ([2} . For 
stellar mass, we use the lLi fc White! (|2009l ) stellar mass func- 
tion at z — 0, and we linearly inte rpolate to higher r edshif t 
stellar mass functions as given in iMarchesini et al.l (|2007|) . 
We consider two different z — halo mass bins: Afhaio ~ 10^^ 
/i"^MQ and lO"'^ /i"^Mq. Using our Af-body simulation, 
we track these halos back in time and determine their masses 
as a function of redshift, assigning them galaxy masses that 
lie on the M«/Mhaio relation at each redshift. Assuming a 
galaxy is quenched at redshift z, we calculate the difference 
between the stellar mass at z and the mean relation at z = 
for each halo mass bin. The earlier a galaxy is quenched, the 
larger the difference with respect to the the typical galaxy 
mass for that halo mass. The currently best-estimated log- 
nor mal scatter in the M<,-Mhaio relation at z ~ is 0.15 — 0.2 
dex (|More et al.ll201Gl : iLeauthaud et all 1201 Ibl ). so the off- 
set in Fig. [14] becomes larger than the scatter for centrals 
quenching prior to z ~ 0.3. 



6 DISCUSSION 

6.1 The Environmental Dependence of Galaxy 
Properties 

In this paper we have demonstrated that the properties of 
galaxies at fixed halo mass are independent of density mea- 
sured on scales larger than the halo radius, and that the 
dependence of the quenched fraction of galaxies on large- 
scale environment can be explained by the variations of the 
halo mass function with environment and by disentangling 
the relative contributions of central and satellite galaxies. 
Our results extends earlier works that demonstrated that 
galaxy properties correlate weakly with large-scale environ- 
ment when small-scale environment is fixed o n a 1 Mpc scale 
(|Kauffmann et al.|[20o3 : JBlanton et al.ll2006l '). However, we 
have argued that, because different galaxies live in different 
halos, it is not possible (or physically motivated) to find a 
single smoothing scale that fully determines galaxy proper- 
ties at all masses. For rich groups and clusters , our r esults 
are in good agreement with lBlanton fc Berlindl (|2007|), who 
found that the the blue fraction of galaxies is independent of 
density at fixed group luminosity, though we do see evidence 
for some change with density for the brightest central galax- 
ies. These results are in good agreement with the results of 
IPeng et al.l (|201H ) , who also find that the environmental de- 
pendence of the red fraction of galaxies is driven primarily by 
satellite galaxies and not by centrals. In detail, the observa- 
tions diverge from this scenario for quenched bright central 
galaxies and in the mean stellar ages of central galaxies on 
the star- forming sequence. But our results support the no- 
tion that halo mass is the dominant quantity governing the 
properties of the galaxies contained within them. 



© 2011 RAS, MNRAS 000.[TH2n 



Halo and Galaxy Formation Histories 13 



6.2 Halo Assembly Bias and the /q-Density 
Relation 

As discussed in §1, the formation liistories of darlc matter 
halos depend on environment. For low-mass halos, old ha- 
los form in high-density environments. For high-mass ha- 
los, old halos form in low-density environments, although 
there is some debate on the amplitude of this correlation. 
If halo and galaxy formation histories are correlated, then 
our measurements are inconsistent with both of these the- 
oretical results. Figs. [7| and [11] clearly demonstrate a lack 
of correlation between /q and p/p for low-mass galaxies. 
Fig. [TJ] demonstrates that the mean stellar ages of quenched 
central galaxies also is constant with p/p. For galaxies be- 
low the knee is the stellar mass function, the ansatz that old 
galaxies live in old halos is not supported by the data. 

However, the quenched fraction of central galaxies in 
Mhost ~ lO'^^'^ h~^ M0 halos show a positive correla- 
tion with environment. This is the mass scale at which 
halos in simulations show little-to-no correlation between 
age and p/p. One possibility is that the environmental de- 
pendence for high-mass central galaxies arises from those 
that have passed through the virial radius of a more mas- 
sive halo (and therefore reside higher density envi ronments) 
but are counted as centrals at the present time (jGill et al.l 
l2005l : lLudlow et al.|[2009l : IWang et al.ll2009h . However, these 
ejected satellites are expected to become less common at in- 
creasing galaxy mass, so the effect is opposite to the trend 
we see here (we will discuss these processes in more detail in 
Papers II and HI). Another possibility is that our assump- 
tion that the central galaxy is the most massive galaxy is 
becoming les s valid at higher galaxy mass, as suggested by 
ISkibba et al.l (|201ll ). and we are confusing high- mass cen- 
trals and satellites. However, in order to masquerade as 
the signal in Fig. IIOI this effect would have to have strong 
environmental dependence, which remains unclear. Addi- 
tionally, at higher halo masses, central and satellite galax- 
ies have increasingly similar quenched fra c tions. Thus, the 
stellar mass dependences of ISkibba et al.l (|201ll ) and cen- 
tral/satellite quenched fractions should cancel out at some 
level. 

Our re s ults have mixed agreement with those of 
IWang et al.l (|2008i ). who used a similar group-finding algo- 
rithm. They found that galaxy groups with redder overall 
colors have higher clustering than groups of the same mass 
but with bluer galaxies, in agreement with our results. How- 
ever, they also found that the effect is stronger in less mas- 
sive groups, opposite to the mass dependence we see here. 
Wang et. al. use color as a proxy for age as opposed to 
(dust-insensitive) Dn4000, though it is unclear why using 
color would lead to a stronger signal. 



6.3 The Stellar to Halo Mass Ratios of Low-Mass 
Galaxies 

Galaxies that are star forming and galaxies that are 
quenched form in halos that growth the same average 
rate. This implies that the stellar mass to halo mass ra- 
tio, Mt/Mhaio will evolve differently for quenched and star- 
forming galaxies. Fig. [T?] demonstrates a Mhaio = lO^'^'^ 
h~^ M0 halo that contains a central galaxy quenched at 
z — 0.6 will have a A/*/Afhaio ratio lower by a factor of 



2.5 at 2 = relative to the mean mean for that halo mass. 
However, weak lensing and satellite kinematical measure- 
ments of the halos masses of low-mass galaxies find that the 
halo masses ar e the same within the e r rors, regardless o f 
color selection JMandelbaum et al.|[2006l : lMore et al.ll201(]| ). 
This is in conflict with the conclusion that galaxy growth 
and halo growth are uncorrelated. There are several possi- 
ble resolutions to this discrepancy. 

(i) The observations of halo masses of low-mass central 
galaxies are biased due the use of color rather than D„4000 
for sample segregation. The low-mass red samples will con- 
tain many intrinsically star-forming galaxies, thus bringing 
the mean halo masses closer together. Lack of significant 
numbers also increase observational errors and low masses. 

(ii) Low-mass galaxies on the red sequence were quenched 
recently, and thus have not had enough time to create a 
signiflcant difference in the M«/Mhaio ratio. There is a clear 
red sequence at z = 1, but the fraction of quenched central 
galaxies at M, < 10 ^"^ is not known. 

(iii) Sub-L» central galaxies do not follow a one-way path 
to the red sequence; galaxies migrate back and forth from 
the blue cloud to the red sequence before permanently be- 
coming quenched. Some blue galaxies at z = may have 
been temporarily quenched previously, while some z = 1 
red galaxies may still grow through star formation later. 
This would minimize the difference between the halo masses 
of red and blue central galaxies at fixed M* . 

Galaxies that are quenched, for example, by AGN 
heating or major galaxy mergers (e.g., Croton et al.l 20061: 
Bower et all l2006l : iHopkins et"ZI l2008l : ISomerville et al] 



2003) could still accrete gas or merge with gas-rich ob- 
jects at some later time, replenishing their fuel supply 
and moving them back into the star-forming sequence. 
In hydrodynamical simulations of galaxy formation, low- 
mass halos ( Mhaio 5, lO'" ^ ^~^ Mp)) accrete their gas in 
a cold state (iKeres et al.l 120051: iDekel fc BirnboimI I2OO6I : 



iKeres et al.ll2009l : iBrooks et al.ll2009f ). Thus previous merg- 
ing or AGN activity that might temporarily quench a galaxy 
will not prevent any new gas from migrating to the center of 
the potential well within a dynamical time, even if the stel- 
lar population is old. However, none of these solutions are 
mutually exclusive, and it is possible that a all three cause 
the measurements presented here. 

6.4 What Quenches Low-mass Central Galaxies? 

The fraction of quenched central galaxies in the deepest 
voids is the same as that in regions around clusters. Thus 
the mechanism by which low-mass central galaxies stop 
their star-fo r matio n must be environmentally independent. 
IWang et al.l (|2009l ) note that many red dwarf galaxies are 
isolated, but propose that the majority of quenched central 
galaxies of this type are ejected satellites, low-mass halos 
that have passed through a group or cluster. The unifor- 
mity of /q°" for low-mass galaxies argues strongly against 
such a model due to the strong environmental dependence 
that would result from such a model. We cannot rule out two 
mechanisms that operate primarily at low and high density, 
respectively, that add to a constant /q°", but Occam's razor 
argues against such a scenario. 

We find that the Sersic indices of these galaxies does 



© 2011 RAS, MNRAS 000.[TH2n 



14 Tinker, Wetzel & Conroy 




-I— I— I— I— p-i— I— I— I— T^ 
SDSS Centrals • 

Halo Growth from; 
z=1.0 

z=0.5 

z=0.3 



f-l 



9.5 10 10.5 

log M,/[h-2Mj 



11 



0.8 



0.6 



X 



B. 0-4 - 



< 



0.2 - 



- 



-| — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



Mh=lx 10'2 



Mh=3x 10" 



_L 



_L 



_L 



_L 



0.2 0.4 0.6 0.8 
quenching redshift 



Figure 12. The ratio of the central quenched fraction between 
high- and low-density environments as a function of stellar mass. 
Filled circles show results from the group catalog, while curves in- 
dicate results using the growth rates of corresponding dark matter 
halos as measured over different redshift baselines. 



Figure 14. The difference in logarithmic stellar mass between 
the 2 = mean A^»-Mhaio relation and the relation for galaxies 
that were quenched at redshift 2, as a function of z. This plot 
assumes that halo growth and galaxy growth are uncorrelatcd, 
thus when a galaxy is quenched the halo it resides in continues 
to grow at the mean rate up to z = 0. Solid and dashed curves 
arc for halos with z = masses of 3 X 10^^ and lO'^^ h~^ Mq, 
which host central galaxies of logAf* = 9.7 and 10.2 h~^M-Q, 
respectively. 



not vary with density, with (n^) ~ 3.5. This implies that 
the mechanism that quenches low-mass centrals also induces 
morphological transformation. It also strengthens the claim 
that the process that halts star formation in centrals at high 
densities also operates at low dens ities. 

Using semi-analytic models, ICroton fc Farrarl (|2008l ) 
find a small population of central galaxies quenched by AGN 
heating, but this process can only effect high-mass halos of 
-Mhaio > 10^^'^ h~^ Mq, roughly an order or magnitude more 
massive than the host halos of the galaxies prob ed here, and 
even m ore massive than the dwarf galaxies in IWang et al.l 
(|2009l ). 

Halo merge rs vary only weakly w ith large-scale dark 
matter density (jFakhouri fc Mall2009l ). implying that the 
galaxy merger rate would also show no environmental de- 
pendence. Major mergers are also likely to cause a transi- 
tion from disk -dominated to sphero idal-dominated morphol- 
ogy (see, e.g.. iHopkins et al.ll2008l ). But it remains unclear 
whether there are enough mergers to fully account for these 
objects. We leave this to a future study (Wetzel et. al., in 
preparation). 

Further analysis of the morphologies and other proper- 
ties of these galaxies, as well as their redshift evolution, will 
shed fight on the mechanisms by which quenching occurs in 
these systems, and whether or not quenching is transient. 



The authors would like to thank Frank van den Bosch for 
supplying the Yang et. al. group catalog, as well as useful 
discussions on the results presented in this paper. 



APPENDIX A: REAL VS. 
DENSITY MEASURES 



REDSHIFT-SPACE 



Here we test the correlation between between large-scale 
galaxy density measured in real and redshift spac e . We use 
the the numerical simulation of I Wetzel fc Whitd (|2010l ) in 
which dark matter (sub)halos are matched with galaxy lumi- 
nosity by the abundance-matching method described in §2, 
and the subhalos ar e populated into a larger volume simula- 
tion as described in [Tinker fc Wetzell (|2010f ). The volume of 
the simulation if 720^ (ft"^ Mpc)^, larger than the Mr < -20 
volume limited sample from Table 1. For each galaxy in the 
simulation, we determine the galaxy density in 10 h~^ Mpc 
spheres in both real space and redshift space, assuming the 
distant observer approximation thus using the 2-axis of the 
simulation as the line of sight. Fig. lAll shows the correla- 
tion and scatter between the two. At this scale, the mean 
redshift-space density faithfully tracks the real-space den- 
sity. More importantly, the scatter is generally < 0.1 dex 
within the extremes of the density field. 



APPENDIX B: GROUP FINDING ALGORITHM 

Here we describe our galaxy group finding algorithm in more 
detail. We define dark matter such that the mean density A 
interior to a radius Ra is 200 times the universal matter 
density p = flm x Pcrit- Here and throughout this paper, Ra 
is in comoving units. Thus the mass and radius of a halo are 
related by 
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Figure 13. Panel (a): D„4000 values for quiescent central galaxies as a function of p/p. The connected dots indicate the mean D„4000 
for different M* bins. The shaded regions indicate the dispersion for each stellar mass bin. Panel (b): The change in the mean quenching 
timescale inferred from (D„4000) in panel (a). Timescales are calculated from stellar population synthesis models described in the text. 
Here we plot At rather than absolute age because it is less sensitive to model assumptions. Connected dots represent the same stellar 
mass bins as in (a). The errorbars indicate the dispersion in At implied by the variance in D„4000. Black curves represent results from 
dark matter halos: the solid black curve is the change in the mean formation time tgs for halos of mass 10^^'^ h^^ Mq (equivalent to 
the halo mass for M* = 10^'^ Mq/H^ central galaxies). Dashed black curve shows At for the oldest 20% of halos at each density bin 
(equivalent to /q for central galaxies with Af* = 10^'"^ Mq/H'^). The dotted curve shows At for halos based on tii2- Panel (c): same as 
(a), but now for star-forming central galaxies (defined as Dn4000 < 1.6). Panel (d): Same as (b) but again for star-forming galaxies. 
In this panel, the dashed curve indicates At for the 80% youngest halos at each density (corresponding to the fraction of star-forming 
central galaxies). Note the change in j/-axis range. 
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Figure Al. Comparison between galaxy density in redshift 
space and real space as measured in a sphere of radius 10 
h~^ Mpc. Results are presented for mock galaxies within a 720^ 
(h^^ Mpc)^^ simulation volume. The solid red curve shows the 
mean relation, and the shaded region shows the dispersion about 
the mean. The solid black curve shows where the two measures 
are equal. 



Ma = ^7r7?iAp. (Bl) 

In comoving units, p is a constant, thus the comoving radius 
of a halo is independent of redshift at fixed mass. Assuming 
the virial theorem, the velocity dispersion of the dark matter 
within a halo is 

2 GMa 



2i?A 



■(1 + ^), 



(B2) 



where G = 4.304 x 10"'' in units of (km/s)^Mpc/M0, and 
the factor of (1 + z) is to convert the comoving radius to 
physics units. 

Dark m atter halos foll o w a 'u niversal' density profile 
described in iNavarro et all |l997i) (hereafter NFW), de- 
fined by a scale radius r^ and a concentration parameter 
ca = Rt^JTs- The 2-D projected density profile of the NFW 
profile is given by 



E(i?) = 2 r, 
with 



5 Pgal f{Rlrs 



/(^) 



i^i^-] 



y X 



if a; < 1 

iix = \ 
if a; > 1 



(B3) 



(B4) 



and 



200 



C200 



3 ln(l + C200) — C2oo/(l + C200)' 

where Pgai is the mean galaxy density in the volume-limited 
sample. Thus, in constructing galaxy groups, we assume that 
halos defined by a mean overdensity with respect to dark 



matter of 200 will also appear at overdensities in the galaxy 
density field with A = 200. 

We assume that the velocity dispersion of satellite 
galaxies follows a Gaussian distribution with dispersion 
given in equation (|B2p . Thus the line-of-sight probability 
of a galaxy being a member of a group is defined by 



p(Ai;) = 



1 



lixa. 



■ exp 



'(A«) 



IgI 



(B6) 



where Ai; = (zgai — 2group)/c and c is the speed of light, and 
the redshift of the group is defined as that of the central 
galaxy. 

All together, the number density contrast of a group 
of galaxies over the background, as a function of both pro- 
jected separation from the central galaxy R, and line-of-sight 
velocity difference from the central galaxy Aw, is expressed 



P{Av,R) 



c T.{R) 

Ho Pgal 



p(Aw) 



(B7) 



For a given galaxy, if P is greater than a background value B, 
then the galaxy is considered to be a member of the group. 
We have calibrated B — 0.5 via our simulation mock cata- 
log to maximize both completeness and purity. We proceed 
through the galaxy catalog in descending order of luminos- 
ity, and if a galaxy is classified as a satellite, it cannot be a 
central galaxy (by definition) and thus cannot have satellites 
of its own. In this algorithm, the brightest (or most massive) 
galaxy in a group is the central galaxy by definition. 

Once galaxies have been assigned as centrals and satel- 
lites in dark matter halos, we recalculate halo mass. But we 
are now interested in deconvolving the galaxy distribution 
into the occupation of galaxies in host halos, we use the 
same abundance matching technique, but now we use the 
total group luminos ity and match tha t with the host halo 
mass function fthe [Tinker et al.l 120081 function). With the 
new halo masses, we go back to the beginning and re-classify 
galaxies as centrals and satellites, iterating the process until 
convergence. 



APPENDIX C: TESTS OF THE GROUP 
FINDING ALGORITHM 

Th e group find e r desc ribed above has been thoroughly tested 
bv lYang et al.l (|200a ). Because we are using a slightly dif- 
ferent implementation of the algorithm, and because we are 
probing density-dependent statistics that have not been ex- 
plored in previous papers, we perform two additional checks 
on our group catalog. 



CI The Halo Occupation Distribution of the 
Group Catalog 

The aim of the group finding algorithm is to associate every 
galaxy in the sample with a dark matter host halo, in the 
same framework as the halo occupation distribution (HOD) . 
For a fixed cosmology, a given HOD maps uniquely onto a 
two-point correlation function. Thus a consistency check on 
the results of the group finder is to compare the observed 
correlation function of galaxies with that predicted by the 
HOD measured from the group sample. 
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Figure Bl. Panel (a): Circles represent the halo occupation of satellite galaxies with M^ < —19 as a function of halo mass as measured 
from the group catalog. The solid curve shows a parametric HOD fitted to these data; the curve that asy mptotes to log Afg^ l = 
represents central galaxies. The dotted line is the power-law fit to N^^i for the same sample of galaxies from lYang et al.ll2008l (their 
fit is restricted to Mijaio i^ 12.3). Panel (b): Points with errors represent measurements of the projected two-point correlation function, 
Wp{rp), from DR7 for A/,. < —19 galaxies. The solid curve is the HOD prediction for Wpivp) using the occupation function from panel 
(a). 



Fig. IBlh shows the the number of sateUite galaxies with 
Mr < —19 as a function of host halo mass. The dotted 
cur ve shows th e pow er-law fit to the same quantity from 
the lYang et al.l (J2008h SDSS DR4 resuhs. They do not in- 
clude results from {Nsat}M < 0.1 host halos in the fit. The 
black curves shows a fit to bo th (NBat)M and (Nc en)M using 
the HOD parameterization in lTinker et al.l (|2008h . Fig.lBTb 
shows the measured projected two-point cor relation function 
o(rp ) from DR7. Here we have used the iLandv fc Szalavl 



(jl993r i estimator for Wp{rp) and the random catalogs pro- 
vided in the VAGC to account for the survey angular selec- 
tion function. The black curve shows the predicted Wp(rp) 
from the HOD. W e use the analytic model described in 
[Tinker et al.l ()2005l ). The two correlation functions are in 
excellent agreement with one another, demonstrating the 
consistency of the halo occupation of the groups. 

C2 Recovery of Galaxy Properties on Mocks 

It is not possible to determine absolute membership of a 
group in redshift space. Group members may have large 
velocities relative to the group center, while non-members 
may be redshifted into the group from coherent infall. Thus, 
some fraction of the group members are not true members 
but are central galaxies that exist within a smaller halo. Al- 
ternately, some galaxies classified as low-mass centrals are 
actually members of a more massive group. (Satellites can 
also be misclassified across groups, but this is comparatively 
rare.) The test presented in Fig. IB II demonstrates that the 
mean number of members is not biased, but because satel- 
lite galaxies have different properties than central galaxies of 
the same magnitude, the mean properties of observed group 



members will be biased. The true quenched fraction of satel- 
lite galaxies is higher than the true quenched fraction of cen- 
trals, thus the misclassification of central/satellite galaxies 
can only have the effect of reducing the observed /q*' and 
increasing /q^". 

To test our group finding algorithm, we apply it to a 
mock galaxy catalog from a 250 h~^ Mpc Af-body simula- 
tion with a particle mass of 10^ h~^ Mp and force resolution 
of 2.5 h'^kpc (see IWhite et al.ll201ol ). which robustly re- 
solves the subhalos that host galaxies we examine here. We 
describe the simulation, subhalo finding, and mock catalog 
generation in detail in Paper HI. The cosmological parame- 
ters of this simulation are nearly identical to those assumed 
in this paper. Galaxy luminosities are assigned to halos and 
subhalos using the subhalo abundance matching technique 
described in §2. As defined in §1, galaxies within subhalos 
are labeled satellite galaxies, while all other galaxies are cen- 
tral. The positions of the galaxies are transformed from real 
space to redshift space using the distant observer approxi- 
mation with the z-axis of the box being the line of sight. 

The left panel in Fig. lCll shows purity and completeness 
of the groups determined via the mock. These statistics are 
broken down into central and satellite galaxies, plotted as 
a function of (intrinsic) host halo mass. Across most halo 
masses, the purity and completeness of central galax;ies is 
over 90%. The purity and completeness of satellite galaxies 
is lower, ~ 80%. The average overall fraction of galaxies 
that are misclassified — labeled as a central when they are 
intrinsically a satellite (or vice versa) — is 10%. 

To test our recovery of the quenched fraction of galax- 
ies, we randomly select 20% of all centrals to be quenched 
and we randomly select 80% of all satellite galaxies to be 
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quenched. These values were chosen to make the difference 
between central and satellite galaxies large in order to maxi- 
mize any bias accrued by the group finder. The impurity and 
incompleteness of the group finding algorithm yields errors 
in the observed quenched fraction of central and satellite 
galaxies. The right panel in Fig. IC II shows the resulting val- 
ues of /q as a function of density, similar to Fig. [7] The 
results of the group finder are shown with the circles with 
error bars. The observed /q'' is reasonably accurate, but 
the observed /q'* is biased low by more than 10%. At higher 
densities, the fraction of satellite galaxies increases, yielding 
a monotonically increasing /q°" with p/p. 

If the fraction of misclassified galaxies is known, then 
this bias can be corrected for statistically. The number of 
quenched satellite galaxies is defined as 



(sat) 
sa-tJQ,obs 



iVsat/, 



^ a^/q '^-t + (1 - 



y)/Q"^ccn, 



(CI) 



where x is the fraction of satellite galaxies correctly identi- 
fied as satellites, and y is the fraction of centrals correctly 
identified as centrals. These two quantities are related by 



(l-a;)Afeat = (l-y)iVccn. 

Equation (|C1|) reduces to 



/, 



(sat) 
Q,obs 



r^fT + {i-x)r^^ 



(C2) 



(C3) 



and the observed quenched fraction for centrals can be sim- 
ilarly written as 



(.(ceil) /H \ ^sat I /.cell 

/q.obs = (1 - yjfci + yjQ ■ 



(C4) 



Solving equation (|C4|) for /q*" and substituting into equa- 
tion (|C3I) . the intrinsic /q** is given by 



/-sat 

/q — 



?//tobs + (^- 



1 ^ J-(cen) 
-^"Q.oba 



xy + {x- l)(l-y) 



(C5) 



Substitution of fq back into equation ([04]) yields the equa- 
tion for /q^". We find that x ~ 0.8 in our A/^-body simulations 
(see Fig. ICll left panel), nearly independent of halo mass 
and large-scale density. To determine the overall quenched 
fraction of centrals, the misclassified fraction of centrals is 
y « O.SiVsat/iVceii = 0.93. To determine the /^™ as a func- 
tion of density, y is calculated at each value of p/p. In prac- 
tice, this procedure overcorrects for /q^" at the densities 
where A^sat > A'cen (roughly p/p > 6, or the rightmost 2 
data points in Fig. IClf) . In this regime, we fix y = 0.85 and 
calculate x through equation (|C2|) . 

The solid curves in right panel of Fig. ICll show the re- 
sults of applying the above procedure to the measured values 

nf f'-"™'' and f''''"' 
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Figure CI. Left panel: The purity and completeness of the group finder vs. halo mass for all galaxies with Mr < —19, determined 
from applying the group finder to a mock galaxy distribution from simulation. Purity is defined by the fraction of centrals/satellites in 
the group finder that were of the same type in the input catalog, while completeness is defined by the fraction of centrals/satellites in 
the input distribution that were assigned to the same type by the group finder. Right panel: Quenched fractions for central and satellite 
galaxies vs. 10 h~^ Mpc overdensity. Dotted lines show the input values, circles show the measurements from the group finder, and solid 
curves show the corrected values from equations l|C5|l and l|C3[l . 
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